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ABSTRACT 

We present LITTLE THINGS (Local Irregulars That Trace Luminosity Ex- 
tremes, The H I Nearby Galaxy Survey) that is aimed at determining what drives 
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star formation in dwarf galaxies. This is a multi-wavelength survey of 37 Dwarf 
Irregular and 4 Blue Compact Dwarf galaxies that is centered around Hl-line 
data obtained with the National Radio Astronomy Observatory (NRAO) Very 
Large Array (VLA). The H l-line data are characterized by high sensitivity (< 1.1 
mJy beam"^ per channel), high spectral resolution (<2.6 km s~^), and high an- 
gular resolution (~6"). The LITTLE THINGS sample contains dwarf galaxies 
that are relatively nearby (<10.3 Mpc; 6" is <300 pc), that were known to con- 
tain atomic hydrogen, the fuel for star formation, and that cover a large range 
in dwarf galactic properties. We describe our VLA data acquisition, calibration, 
and mapping procedures, as well as Hi map characteristics, and show channel 
maps, moment maps, velocity-flux profiles, and surface gas density profiles. In 
addition to the H I data we have GALEX UV and ground-based UBV and Ha 
images for most of the galaxies, and JHK images for some. Spitzer mid-IR im- 
ages are available for many of the galaxies as well. These data sets are available 
on-line. 

Subject headings: galaxies: irregular — galaxies: star formation — galaxies: ISM 
— galaxies: kinematics and dynamics — galaxies:structure 

1. Introduction 

Dwarf Irregular (dirr) galaxies are the most common type of galaxy in the local universe 
(see, for example. Dale et al. 2009). They are also the closest analogs to the low mass 
dark matter halos that formed after the Big Bang. In the ACDM model, it is in these 
entities that the first stars formed. Yet, we do not understand what drives star formation 
on galactic scales even in nearby dwarfs, the simplest, most pristine local environments (see, 
for example. Hunter 2008, Bigiel et al. 2010). To remedy this situation we have assembled 
multi-wavelength data on a large sample of relatively normal, nearby gas-rich dwarf galaxies, 
tracing their stellar populations, gas content, dynamics, and star formation indicators. This 
project is called LITTLE THINGS (Local Irregulars That Trace Luminosity Extremes, The 
Hi Nearby Galaxy Survey), and it builds on the THINGS project, whose emphasis was on 
nearby spirals (Walter et al. 2008). 

LITTLE THING^was granted ~376 hours with the Very Large Array (VLA0) in B, 



^^http: / / www.lowell.edu/users / dah/littlethings /index. html 

1^ The VLA is a facility of the National Radio Astronomy Observatory (NRAO). The NRAO is a facility 
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C, and D configurations to obtain H l-line data of 21 dirr and Blue Compact Dwarf (BCD) 
galaxies. We added another 20 galaxies from the VLA archives for a total sample of 37 dIrr 
and 4 BCD galaxies covering a large range of dwarf galactic parameter space. Our VLA data 
emphasize deep, high spatial (6") and spectral (<2.6 km s^^) resolution in order to reveal 
clouds, shells, and turbulence in the interstellar medium (ISM) that could be important for 
creating regions of higher density star-forming gas. We are combining the H l-line data with 
ultraviolet (UV), optical, and infrared (IR) data to address the following questions: 

What regulates star formation in small, gas-rich galaxies? Studies of the star formation 
laws in galaxies show that star formation is very inefficient in dwarf galaxies and in the 
outer disks of spirals (Bigiel et al. 2010; Ficut-Vicas et al., in preparation). Furthermore, 
the star formation efficiency decreases as the gas density decreases. How then do atomic 
hydrogen clouds form in the low gas density environments of dwarfs that ultimately lead to 
star-forming molecular clouds? 

What is the relative importance of sequential triggering by previous generations of stars 
for star formation in dwarf galaxies? One generation of stars can trigger the formation of the 
next by rearranging the gas through winds and supernova explosions (see, for example, Opik 
1953, Gerola et al. 1980, Comins 1983, Dopita et al. 1985, Mori et al.l997, van Dyk et al. 

1998), but how important is this process? Observations of dwarfs show a better correlation 
between the star formation rate and the \/-band surface brightness, which emphasizes Gyr 
old stars, than with any other measure, including the average gas density (Hunter et al. 
1998). Does this relationship imply a significant causal relationship in which existing stars 
are important for triggering new star formation? 

What is the relative importance of random turbulent compression cloud formation in 
dwarf galaxies ? Turbulence can account for a wide range of phenomena that are indirectly 
related to star formation: the power law luminosity functions of Hll regions, the power 
law mass functions of clouds and clusters (Mac Low & Klessen 2004), the shape of the 
power spectra of Hi (Stanimirovic et al. 1999, Elmegreen et al. 2001, Dickey et al. 2001, 
Zhang et al. 2012a) and Ha (Willett et al. 2005), the log-normal distribution of the Ha 
intensity in individual pixels in Ha images of dwarf galaxies, (Hunter & Elmegreen 2004), 
the correlation between region size and the star formation time scale (Efremov & Elmegreen 
1998), and hierarchical structures in the ISM (Bastian et al. 2007). But how important is 



of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
These data were taken during the upgrade of the VLA to the Expanded VLA or EVLA. In this paper we refer 
to the instrument as the VLA, the retrofitted antennae as EVLA antennae, and non-retrofitted antennae as 
VLA antennae. This emphasizes the hybrid nature of the instrument and distinguishes it from the far more 
powerful Jansky VLA or JVLA it has become since 2012. 
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turbulence in the formation of star-forming clouds, particularly in the realm of sub-critical 
gas density? 

What happens to the star formation process in the outer parts of disks? We find that 
many dwarfs have exponential stellar surface brightness profiles with a break and a change 
in slope. These breaks are similar to what are seen in the outer parts of spirals (e.g., de 
Grijs et al. 2001; Kregel et al. 2002; Pohlen et al. 2002; Kregel & van der Kruit 2004; Pohlen 
& Trujillo 2006; Erwin et al. 2008; Gutierrez et al. 2011) and disks at high rcdshift (Perez 
2004). The break could imply a change in the star formation process at the break radius, 
but what is really happening there? 

Here we describe the LITTLE THINGS project and present the H I data. Science papers 
that are part of this project include: 1) "Velocity Dispersions and Star Formation Rates in 
the LITTLE THINGS Dwarf Galaxies (Cigan et al., in preparation), 2) "Star Formation 
Laws in LITTLE THINGS Dwarfs: The case of DDO 133 and DDO 168" (Ficut-Vicas et 
al. , in preparation) , 3) "Deep Radio Continuum Imaging of the Dwarf Irregular Galaxy IC 
10: Tracing Star Formation and Magnetic Fields" (Heesen et al. 2011), 4) "Deep Radio 
Continuum Observations of the Nearby Dwarf Irregular Galaxy IC 10. I. Nature of the 
Radio Continuum Emission" (Heesen et al., in preparation), 5) "Surface Brightness Profiles 
of Dwarf Galaxies. I. Profiles and Statistics" (Herrmann et al. 2012), 6) "Surface Bright- 
ness Profiles of Dwarf Galaxies. II. Color Trends and Mass Profiles" (Herrmann et al., in 
preparation), 7) "The Stellar and Gas Kinematics of the LITTLE THINGS Dwarf Irregu- 
lar Galaxy NGC 1569" (Johnson et al. 2012), 8) "The Stellar Kinematics of the LITTLE 
THINGS Dwarf Irregular Galaxies DDO 46 and DDO 168" (Johnson et al., in preparation), 
9) "Deep 6 cm Radio Continuum Observations of the Nearby Dwarf Irregular Galaxy IC 
1613" (Kitchener et al., in preparation), 10) "High-resolution Mass Models of Dwarf Galax- 
ies from LITTLE THINGS" (Oh et al, in preparation), 11) "Star Formation in IC 1613" 
(Simpson et al., in preparation), 12) "Outside-in Shrinking of the Star-forming Disk of Dwarf 
Irregular Galaxies" (Zhang et al. 2012b), and 13) "Hi Power Spectrum Analysis of Dwarf 
Irregular Galaxies" (Zhang et al. 2012a). Other science that is in progress includes, but 
is not limited to, studies of the following: 1) Star Formation Processes in Blue Compact 
Dwarfs - Ashley, 2) Molecular Cloud Structure and Dust at Low Metallicity - Cigan, 3) Star 
Formation Laws of the full sample of galaxies - Ficut-Vicas, 4) Surface Brightness Profiles 
and their relationship to Hi gas and two dimensional images - Herrmann, 5) The Role of 
Turbulence in Star Formation in Dwarfs - Hunter, 6) A Deep 6 cm Radio Continuum Survey 
of LITTLE THINGS - Kitchener, and 7) Porosity of the Interstellar Medium in Dwarfs - 
Simpson. 
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2. The Survey 

2.1. Sample selection 

The LITTLE THINGS sample of 37 dirr galaxies and 4 BCDs was drawn from a larger 
multi- wavelength survey of 94 relatively nearby dIrr galaxies and 24 BCD systems, as well as 
18 Sm galaxies (Hunter & Elmegreen 2004, 2006). A large sample is necessary to disentangle 
complex processes in surprisingly complicated galaxies. This larger sample covers nearly the 
full range of galactic parameters seen in dwarf galaxies including luminosity (9 magnitudes in 
absolute magnitude My), current star formation rate normalized to the area within one disk 
scale length (a range of a factor of 10^), relative gas content (a factor of 110 in Mhi/Lb), and 
central surface brightness (7 magnitudes in fi^, the central \^-band surface brightness). The 
sample does not include obviously interacting galaxies because we are investigating internal 
processes, and it only includes dwarfs with substantial amounts of atomic hydrogen so that 
they could be forming stars. 

For LITTLE THINGS, we excluded galaxies with distances >10 Mpc in order to achieve 
a reasonably small spatial resolution in the VLA Hi maps (<300 pc at 10 Mpc in B array), 
although the distance to our most distant galaxy was revised later to 10.3 Mpc. The average 
distance to the galaxies in our sample is 3.7 Mpc. In addition we excluded galaxies with 
W20 >160 km s~^ in order to fit the galaxy emission comfortably in the bandwidth available 
for our desired velocity resolution, leaving room for determining continuum emission. The 
parameter 14^20 is the full-width at 20% of the peak of an integrated H I flux- velocity profile 
(Bottinelli et al. 1990). One galaxy originally in our sample, NGC 1156, was eliminated 
when we discovered that the archival data contained no line-free channels. We also excluded 
galaxies with high inclinations (>70°) so that we can better resolve structures in the disks. 
We then insisted that the LITTLE THINGS sample cover the range of integrated properties 
seen in the full optical survey, especially My and star formation rate, including galaxies with- 
out any measurable Ha emission. A comparison of the properties of the LITTLE THINGS 
sample with the larger parent sample is shown in Figure [T| Our sample includes 19 dIrr 
galaxies and one BCD galaxy in the VLA archives having B/C/D-array observations at 2.6 
km s~^ channel separation or better. Our new VLA observations are of 18 dIrr galaxies and 
3 BCDs in the B/C/D arrays. In total, 19 galaxies have Hl-line data at 1.3 km s~^ channel 
separation. 

One of our 41 galaxies, NGC 6822, which we observed at 1.3 km s~^ channel separation, 
is highly extended compared to the primary beam. Those data require special treatment 
and calibration and mapping are in progress. Below, we list the data we obtained on NGC 
6822, but we do not include it in the discussions of the H I maps and cubes. 
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Table [T] lists the sample galaxies and several basic properties, including other names by 
which the galaxies are known, distances. My, a measure of size, foreground reddening, star 
formation rates, and oxygen abundance. Here we elaborate on a few of the quantities: 1) 
Distance. In choosing a distance for each galaxy, we considered the references for distance 
determinations given by the NASA/IPAC Extragalactic Database (NED). In order of prefer- 
ence, we chose distances determined from Cepheids or RR Lyrae stars and from the tip of the 
red giant branch. The references for the distance used for each galaxy are given in Table [TJ 
If none of those determinations was available, we used the Hubble flow distance determined 
from the galaxy's radial velocity corrected for infall to the Virgo Cluster (Mould et al. 2000) 
and a Hubble constant of 73 km s~^ Mpc"^. 2) Size measurement. For a measure of the size 
of the galaxy, we used the Holmberg radius, Rh, which was originally measured by Holmberg 
(1958) at a photographic surface brightness of 26.5 mag in 1 arcsec^. Here we have used the 
Rh measured by Hunter & Elmegreen (2006), who converted from photographic magnitudes 
to 5-band magnitudes using the recipe of de Vaucouleurs et al. (1976). For our galaxies this 
corresponds to a i?-band surface brightness of between 26.60 and 26.72 mag in 1 arcsec^ with 
an average of 26.64±0.002 mag in 1 arcsec^. 3) Star formation rate (SFR). Star formation 
rates are taken from Hunter & Elmegreen (2004). These were determined from integrated Ha 
luminosities and are normalized to the area encompassed by one disk scale length. We use 
this characteristic area to normalize the SFRs in order to compare galaxies of different size, 
but this is not the area that encompasses Ha emission. Some galaxies have Ha emission that 
extends less than and others more than one disk scale length. 4) Oxygen abundances. As a 
measure of the metallicity of the galaxy we list oxygen abundances, 12-|-log 0/H, taken from 
the literature and determined from H ll regions. For galaxies without H ll region abundance 
measurements, we used Mb and the luminosity-metallicity relationship of Richer & McCall 
(1995). The abundance determined this way, listed in parentheses, is very uncertain, but it 
gives an indication of the expected relative metallicity of the galaxy. The references for the 
abundances are given in the Table. For comparison the solar oxygen abundance value is 8.69 
(Asplund et al. 2009). 



2.2. VLA observational philosophy 

Our new VLA observations are characterized by deep observations with high angular 
and high spectral resolution. To achieve the angular resolution and sensitivity we needed, 
we used VLA B/C/D-array configurations with 12/6/2 hr integration times per galaxy. The 
combined B/C/D arrays sample the galaxies at 6" (110 pc at 3.7 Mpc, the average distance 
of our sample; for ROBUST- weighting). This resolution shows clouds, shells, and turbulent 
structures that are important for star formation. The B array, which is the most extended 
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configuration of the three we used, requires the most time because of its reduced surface 
brightness sensitivity, but the high resolution it provides is essential. Typical clouds and 
holes are of order 100 pc in diameter, and they would be lost at the resolution of the data 
obtained in the more compact C and D arrays (the resolution of C+D is ~20", or 360 pc 
at 3.7 Mpc). However, the C and D arrays are important because they reveal the extended, 
low-density gas around clouds, shells, and star formation. Low spatial resolution maps also 
trace the low density H I far beyond the stellar disks. The D array can detect structure up 
to a scale of 15'. 

For good velocity resolution, wc chose a channel separation of 1.3 km s~^ and 2.6 km 
s^^, both for new observations and for data wc drew from the archives. The choice of channel 
separation was dictated by the need to contain galaxy emission within the bandwidth, leaving 
line-free channels on either side for continuum subtraction. Twenty-two galaxies have only 
Hanning-smoothed data cubes with 2.6 km channel separation, and 19 have small enough 
velocity ranges to be observed and imaged with 1.3 km s~^ channel separation. These 
relatively small channel separations enable the resolution of kinematic structures in the 
galaxies, especially the deconvolution of warm and cool components with velocity FWHM 
as small as 7-9 km s~^ (Young et al. 2003). 

In our sample of galaxies requiring new VLA observations, 13 had no previous VLA 
observations, and 8 had some data in the archives. Altogether, there were 21 galaxies that 
needed B array observations, 15 that needed C, and 17 needed D configurations. Seven 
of our galaxies arc part of the VLA-ANGST (ACS [Advanced Camera for Surveys] Nearby 
Galaxy Survey Treasury) large project (Ott et al. 2010; project ID A0215) that began at 
the same time as LITTLE THINGS: CVnIdwA, DDO 70, DDO 75, DDO 155, DDO 187, 
NGC 4163, UGC 8508. The VLA-ANGST integration times were generally shallower than 
those of LITTLE THINGS, but these galaxies were observed as part of VLA-ANGST rather 
than LITTLE THINGS and we use those data here. In addition eight of our galaxies taken 
entirely from the archives were part of the THINGS (Walter et al. 2008) sample (DDO 
50, DDO 53, DDO 63, DDO 154, M81dwA, NGC 1569, NGG 2366, NGG 4214). For the 
THINGS galaxies, we retrieved their uv-calibrated data. 

There are two galaxies that arc exceptions to our combination of the B/C/D array 
configurations. These galaxies are DDO 216 and SagDIG. The data for these galaxies were 
taken from the archives where only C and D array configuration data were available. 
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2.3. Ancillary data 

We also have a rich multi- wavelength collection of ancillary data on the LITTLE THINGS 
galaxies. Our data sets are summarized in Table [2j We have two tracers of the products 
of recent star formation. First, Ha images obtained at Lowell Observatory reveal the pres- 
ence of massive stars formed over the past 10 Myrs through the ionized nebulae of their 
natal clouds (Hunter & Elmegreen 2004). Second, ultraviolet images at 1516 A and 2267 
A obtained with the Galaxy Evolution Explorer satellite {GALEX; Martin et al. 2005) show 
the bright ultraviolet starlight from stars formed over the past 100-200 Myrs (Hunter et al. 
2010). The ultraviolet images are a more reliable tracer of recent star formation especially 
in the outer disks of galaxies where star formation is more sporadic and lower gas densi- 
ties make Hll regions harder to detect. The GALEX data were processed with the GR4/5 
pipeline. 

Older stars are traced with UBV images obtained at Lowell Observatory (Hunter & 
Elmegreen 2006). For on-going star formation, these are sensitive to stars formed over the 
past 1 Gyr. For 17 of the galaxies, we also have some combination of JHK images obtained 
at Lowell Observatory, and for 36 galaxies there are Spitzer3.Q fim images that we have pulled 
from the archives. The near and mid-IR images in particular highlight the star formation 
activity integrated over a galaxy's lifetime. 

Dust plays an important, although poorly understood, role in the star formation process. 
FIR continuum emission traces this part of the ISM. Spitzer 5.8 fim and 8.0 /xm images 
trace polycyclic aromatic hydrocarbon (PAH) emission (plus hot dust continuum emission). 
Spitzer MIPS images of 35 of the galaxies, obtained as part of the Spitzer Infrared Nearby 
Galaxies Survey (SINGS), Local Volume Legacy (LVL), and other programs, can in principle 
give the dust continuum emission and distribution at 24 /im, 70 /im, and 160 fim, although 
exposures tend to be shallow. We are also mapping a subset of galaxies in the FIR with 
bolometer arrays on the Atacama Pathfinder Experiment (APEX) sub-millimeter telescope 
in collaboration with M. Rubio and with Herschel (Pilbratt et al. 2010; an 0T2, priority 1 
program) in collaboration with S. Madden. These data are important for understanding the 
heating balance in the ISM. 

Molecular clouds are a key component of the ISM and the component most directly 
connected with star formation. We are exploring this part of the star formation process 
through case studies. This has included Combined Array for Research in Millimeter-wave 
Astronomy (CARMA) observations of CO line emission in NGC 3738, a dwarf with intense 
star formation in the center of the galaxy. With APEX we mapped CO in WLM, a typ- 
ical star-forming dwarf, in collaboration with M. Rubio. We are also using Herschel (an 
OTl, priority 1 program) to map photo-dissociation regions (PDRs) in five of the LITTLE 
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THINGS galaxies in collaboration with M. Kaufman and S. Madden. These data will allow 
us to characterize molecular clouds at low metallicity. 

3. VLA observations and data reduction 
3.1. Observing 

The LITTLE THINGS (project AH927) observations began with the B-array 2007 De- 
cember 19 and, for B-array, ended 2008 February 10. The observations were done in two 
6-hour segments for each galaxy. Some observations were repeated when technical hiccups 
occurred. C- array observations followed in the spring of 2008, and D-array observations in 
the summer of 2008. All of our new observations were in hand by the beginning of Au- 
gust 2008. The observations are listed in Table [3j Each array configuration with which a 
galaxy was observed is listed separately. Multiple observing sessions for each configuration 
are identified by the dates on which the observing sessions began. 

Our new observations used a bandwidth of 0.8 or 1.6 MHz (nominally 150 and 300 km 
s~^ wide), depending on the known velocity range of emission from each galaxy. We used 
two intermediate frequencies (IFs; 2AC) and 256 channels. 

Our new data were not Hanning smoothed when the observations were made. Thus, 
the galaxies for which we were targeting a spectral resolution of 2.6 km s~^ (12.21 kHz) were 
observed with a channel separation of 1.3 km s~^ (the spectral resolution is 1.4 times this), 
and in the reduction stage the data were Hanning-smoothed to obtain a channel separation 
and resolution of 2.6 km s~^. Similarly, galaxies for which we targeted a spectral resolution of 
1.3 km s~^ (6.10 kHz) were observed with a channel separation of 0.65 km s^^ (3.05 kHz) and 
later Hanning-smoothed to obtain a channel separation and resolution of 1.3 km s~^. Most 
of the data from the archives were on-line Hanning smoothed. The exception is SagDIG, 
which was observed with a channel separation of 3.05 kHz and no Hanning-smoothing. For 
that galaxy, we produced a cube with the native resolution (1.4 x the channel separation of 
0.65 km s~^), and another cube with Hanning smoothing to produce a channel separation 
and resolution of 1.3 km s~^. 

Because of the transition of the VLA to EVLA during the course of the survey, sensi- 
tivity changes, phase drifts, and phase jumps were expected. To deal with this required that 
careful calibration be built into the observations. Our observations in general began and 
ended with a primary flux calibrator, usually 3C286, allowing extra time at the beginning 
in case of a cable wrap problem in moving to our galaxy position. It was also recommended 
to add a short dummy scan to make sure the antennas were properly initialized at the be- 
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ginning of the observing session. The primary cahbrator was used to set the flux scale and 
for determining the response function across the bandpass. In between these observations 
we alternated between the secondary calibrator and the galaxy with approximately 12-25 
minutes on source at any one time. The shorter on-source dwell times were in the shorter 
observing sessions. Because of our narrow bandwidth, when offered with a choice, we pre- 
ferred a strong calibrator, even if it was at a radial distance that was larger than that of a 
more nearby, but weaker one. This was possible because of the phase stability of the VLA 
at 20 cm even in the B-configuration. Even so, the secondary calibrators were always within 
about 12° of a target galaxy. The secondary calibrator was used to determine the variations 
of phase and amplitude with time during the observing session. The time spent on source 
in each configuration is given in Table [3] 



3.1.1. The Aliasing Problem 

During the course of our observations, the VLA was in transition to the EVLA, and 
antennas were being outfitted with new receivers and electronics one by one. At the beginning 
of our observations 46% of the antennas in the array had been converted, and by the end 62% 
of the array consisted of EVLA antennas. To make the digital signal from EVLA antennas 
compatible with the old analog VLA correlator, the signal on the new receivers was passed 
through an electronic filter. However, this conversion caused power from sidebands to be 



aliased into the lower 700 kHz of the observing bandwidth (the "aliasing problem" ^^). In 
spectral line mode, the line emission would be unaffected as long as the bandwidth used was 
at least 1.5 MHz and the spectral line observed was narrow enough so that it fell beyond the 
affected 0.7 MHz or the line could be shifted away from the band center. However, even so, 
continuum emission was aliased, with phase information scrambled making it resemble noise, 
the effect being strongest in channels closest to the baseband edge. This is most noticeable on 
EVLA- EVLA baselines, and to a lesser extent on EVLA- VLA baselines; VLA- VLA baselines 
were not affected (see Figure |2]). 



This problem had a serious impact on LITTLE THINGS, which used narrow band- 
widths. A key element of our project is high spectral resolution imaging at high sensitivity 
of these faint, low rotation velocity objects. Thus, our galaxies were being observed with 
small bandwidths (0.78 MHz or 1.56 MHz) in order to achieve high spectral resolution (1.3 
km s~^ or 2.6 km s~^). We examined our data to understand the consequences of the 
aliasing problem and how we might deal with it. We took 1.56 MHz bandwidth data of 



^^http://www. via. nrao.edu/astro/guides/evlareturn/aliasing.shtml/ 
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one of our galaxies and fully calibrated and mapped it, running extensive tests. We clearly 
saw increased noise from the lower band edge to the middle of the passband when EVLA- 
EVLA and VLA-EVLA baselines were included (Figure |2]). In the end, we were able to 
satisfactorily calibrate the data only by eliminating the EVLA-EVLA baselines. Subtracting 
the continuum from the line data using a fourth order polynomial, as had originally been 
suggested, resulted in unacceptable amplitude variations across the bandpass. An alternative 
scheme using unaffected baselines to gauge the level of the continuum worked at some level, 
but continuum remained in the map at unpredictable levels. Furthermore, the noise added 
to the overall dataset by the EVLA-EVLA baselines was at such an extent that the noise in 
a data set containing all baselines and the noise in the same data set without EVLA-EVLA 
baselines are at virtually the same level. This means that, even if it were possible to calibrate 
properly the EVLA-EVLA baselines, this would not reduce the noise level. 

Because the EVLA-EVLA baselines constituted 22% of the total in B-array, we requested 
a 22% increase in our B-array observing time, or 55.5 hours, in order to make up for the lost 
EVLA-EVLA baselines and achieve our original target signal-to-noise. This was granted in 
the form of dynamically scheduled time, in blocks of 1-2 observing sessions (2.5-3 hr) per 
galaxy. Similarly, time for C and D array observations were increased in proportion to the 
percentage of EVLA-EVLA baselines lost, and this extra time was dynamically scheduled 
as well. A consequence of the dynamic scheduling was that some of the make-up time 
came during array configuration transitions and not when the instrument was in a "pure" 
configuration (see Table [s]). 

We were worried that the loss of the EVLA-EVLA baselines might change the uv- 
coverage that we expected. Comparing the uv-coverage with and without EVLA-EVLA 
baselines removed, however, showed that the antennas were being retrofitted in a pattern 
that affected the uv-coverage only slightly and caused sidelobe structure that could be dealt 
with in the imaging stage by employing a suitable CLEAN algorithm (see Appendix B). 

3.1.2. Doppler Tracking 

Another consequence of the VLA-EVLA transition was that Doppler tracking was not 
possible. Because the LO chains for EVLA and VLA front-ends are different, frequency 
changes can cause phase jumps between EVLA and VLA antennas. Therefore, observing 
at a fixed frequency was recommended. This meant that the frequency of the observations 
changed by small amounts during the course of an observing session. We chose our central 
frequency for the expected date and time of the observations, and checked to ensure that 
the galaxy would not move out of the channel during the observations. Trickier were the 
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dynamically scheduled observations. For those, we chose a date in the middle of the range 
of possible observation dates. In no case did the galaxy emission move out of the channels, 
but this added an additional concern in setting up the observe files and an additional step 
in mapping the data. 

3.1.3. Avoiding Milky Way Emission 

Another problem arose for six galaxies with radial velocities that are similar to those 
of Hi in the Milky Way in the direction of the galaxy, potentially leading to Milky Way 
emission in the passband. These galaxies include DDO 165, DDO 167, DDO 187, NGC 
4163, NGC 6822, and UGC 8508. (NGC 1569 was observed before the transition to EVLA 
as part of THINGS.) The standard approach in these cases is to observe both primary and 
secondary calibrators at velocities offset from that of Galactic Hi in their direction. This 
is accomplished by applying a frequency offset to a higher and lower central frequency, 
symmetrically with respect to the originally intended central frequency. The offset needs 
to be large enough to shift the entire bandwidth beyond the range where Galactic Hi is 
present. However, because phase jumps were expected to accompany frequency shifts with 
the EVLA antennae, we observed only the primary calibrator with a positive and negative 
offset, ensuring it was unaffected by Galactic foreground Hi. Because no absolute phase 
calibration is required for the flux and bandpass calibration, the phase jumps accompanying 
any frequency shifts were irrelevant for the primary calibrator. For NGC 6822, there were 
two flux calibrators — 3C286 and 3C48, and both were offset in frequency. The offset was 3 
MHz for 4 of the galaxies, and 3.5 MHz for DDO 187 and 2.5 MHz for NGC 6822. 

Unfortunately, in calibrating these data we found a problem with some of the calibrator 
observations on one side of the central frequency. From discussions with Ken Sowinski at 
NRAO we realized that this was due to reassignment of the Local Oscillator (LO) settings 
by the on-line system that tunes the receivers, which meant the frequency range of the 
incoming signal fell partly beyond the range selected by the front end fllter. We, therefore, 
had to throw out the corrupted observation of the calibrator and use only the frequency 
offset observation on the side which was not corrupted. This affects the flux calibration at 
a level of at most 0.1% and can therefore be safely ignored. 
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3.2. Calibration 

The VLA data were calibrated and mapped by many members of the team. Usually 
one person would calibrate all of the data sets for a galaxy, and similarly one person would 
usually complete the mapping of a galaxy. As data sets were calibrated or the ensemble of 
data sets for a galaxy was mapped, the data were uploaded to an internet-accessible hard disk 
provided and maintained by NRAO, making it available to the rest of the team and the next 
person who would work on the data. The final map cubes were then examined by designated 
"checkers" who inspected all of the data, looking for signs of any problems or inconsistencies. 
A galaxy was not considered "done" until it had passed a thorough inspection, and any and 
all problems were resolved. 

In order to ensure uniform treatment of the data, we created "recipes" for data cal- 
ibration and mapping that were carefully followed. Since the THINGS survey (Walter et 
al. 2008) had already been through this process and two of our team members were part of 
THINGS, we started with their recipe and updated it to reflect changes in software, com- 
plications due to the aliasing problem, and problems resulting from the frequency shifts for 
observations that could be contaminated by Galactic Hi. In order to ensure a homogeneous 
data set, we treated the archival data the same way that we treated the new data. All of 
the reductions were done in the Astronomical Image Processing System (AIPS) developed 
by NRAO, version 31DEC08. We used the flux scale of Baars et al. (1977). 

We outline the salient calibration steps in Appendix A, and the recipe itself is available 
on the public web site discussed in ^ After every action on the data, we had at least 
one, and usually multiple, quality control checks. We outline the salient calibration steps, 
including those to check the quality of the calibration, in Appendix A. The recipe itself is 
available on the public web site discussed in §5. For the cases involving multiple frequencies, 
most steps had to be run for each frequency. 

A few circumstances yielded special calibration challenges: 1) Strong winds and many 
technical problems resulting in phase jumps meant that considerable data had to be elimi- 
nated, limiting the data obtained 2008 February 4, affecting DDO 101 (B array). 2) There 
was a clock error on data taken in B array 2008 January 5 and 7 and only VLA- VLA base- 
lines could be used. The galaxies that were affected were DDO 46, DDO 52, and DDO 133. 
A similar problem on 2008 January 24 rendered all data from EVLA antennas unusable, 
affecting DDO 87 and DDO 167. 3) All of the EVLA antennas had to be flagged during 
the archive 2006 B array observations of VII Zw 403. Four antennas were involved. 4) Solar 
interference affected the observations in compact configurations of some galaxies, usually 
those data <1 kA: DDO 43 (C and D arrays), DDO 47 (D array), DDO 52 (D array), DDO 
87 (C array), F564-V3 (D array), IC 1613 (C and D arrays), NGC 3738 (D array), NGC 
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4163 (D array), WLM (B and C arrays), and VII Zw 403 (C and D arrays). 5) There were 
problems with the phases in the right polarization (RR), and some of the left (LL), of the 
primary phase calibrator of DDO 43 (C array) in the archival data obtained in 2000 April. 
To deal with this, we performed the bandpass calibration on the secondary instead of the 
primary calibrator, flagged most of the left polarization data on the primary calibrator, kept 
only the left polarization in creating the new Channel (an average of the inner 75% of the 
frequency channels), ran the usual calibration on the primary and secondary calibrators to 
get the flux of the secondary calibrator, ran the cahbration CALIB only on the secondary, 
and applied the calibration using CLCAL to the new Channel data. 6) There was no model 
for the primary calibrator 3C295 used with VII Zw 403 in D array in the archival 1997 
data. However, the calibrator is not resolved in D array, so we used our observations of the 
calibrator instead of a model. Similarly, the secondary calibrators for archival observations 
of IC 1613 and of WLM in B array were spatially resolved. For IC 1613, we used the C 
array image of the calibrator as a model for calibrating the B array with a maximum limit 
on the uv range of 15 kA. For WLM, the effect was small enough to ignore. 7) One galaxy, 
DDO 101, was self-calibrated in C and D arrays using a 2.3 Jy continuum source in the 
held. 8) The secondary cahbrator, 0925 + 003, of the archive data on DDO 70 had poor 
signal-to- noise (flux of ~0.6 Jy). 9) For IC 10, the D array observation that we used is the 
central pointing (qual=13) of a 25-pointing mosaic (Wilcots & Miller 1998). 



3.3. Imaging 

The mapping recipe picks up where the calibration leaves off, beginning with extracting 
the calibrated uv-data for just the galaxy from the uu-data set (split). The individual data 
sets taken in B, C, and D-configuration, were combined into one large uv-file, ready to be 
Fourier transformed to create maps in the image plane containing structure across a range of 
angular scales, from a few arcseconds up to several arcminutes. The spectra from different 
data sets were aligned in frequency space (cvel), and the uv-daia sets combined (dbcon, 
with no adjustments to the relative weights). We subtracted continuum emission in the 
uv-p\ane (uvlsf), followed by mapping. 

The commonly used mapping/CLEAN algorithms for VLA data approximate all emis- 
sion structure as the sum of point sources. This is clearly an unrealistic approximation for 
extended galactic emission. In addition, the resulting map is the sum of the clean compo- 
nents that are made with the clean synthesized beam and the residuals that are a product 
of the dirty beam. So, for mapping, we were excited by the advantages that multi-scale or 
multi-resolution CLEAN (MS-CLEAN) might bring. This concept, developed by Wakker & 
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Schwarz (1988), solves simultaneously for point sources and Gaussians with pre-determined 
sizes (Cornwell 2008). This has been implemented in AIPS under IMAGR by Eric Greisen 
(Greisen et al. 2009). However, there has been limited experience with using this algorithm 
in the community (Rich et al. 2008). Therefore, before we could adopt it as our standard, 
we invested considerable time experimenting with the algorithm in order to understand the 
consequences to the data, in particular comparing MS-CLEAN maps against those produced 
with standard CLEANing methods (for more details see Appendix B). We concluded that 
the MS-CLEAN algorithm gave superior maps and that we understood it well enough to 
adopt it for making the LITTLE THINGS map cubes. We then finahzed our mapping recipe 
around the AIPS MS-CLEAN option in iMAGR (see details in Appendix B). 

We produced both a Natural-weighted image cube, which we expected would bring out 
extended lower surface brightness emission at the expense of angular resolution and beam 
profile, and one weighted with a ROBUST value of 0.5 (Briggs 1995), which produces a better 
behaved synthesized beam at a resolution that is close to that obtained with Uniform- 
weighting. Usually, the ROBUST- weighted cube was cleaned to 2xrms and the Natural- 
weighted cube to 2.5xrms. These levels were chosen after extensive tests that showed that 
the flux that was recovered by MS-CLEAN was not a strong function of CLEAN level, as long 
as that level was within factors of a few of the optimum level. Cleaning deeper than 2xrms 
gave noticeable CLEAN artifacts; cleaning down to 3-4xrms left significant residual fiux in 
the different scale maps. We converged on a level of 2xrms. However, the Natural- weighted 
cubes cleaned to 2xrms showed signs of slight over-cleaning ("holes" dug in the background). 
For the Natural- weighted cubes we found that 2.5xrms was the best compromise. Each map 
cube was carefully inspected for signs of over- or under-cleaning. 

The maps of two galaxies, Mrk 178 and VII Zw 403, required small tweaks to this 
scheme. For Mrk 178, our standard CLEAN to 2xrms with the ROBUST- weighting produced 
a field at the 15" scale (see Appendix B) that appeared to be over-cleaned. We, therefore, 
cleaned that scale to 2.5xrms, but the other scales to 2xrms. For VII Zw 403, a CLEAN 
depth of 2xrms resulted in being able to see the CLEAN box in the first two fields (0", 15"), 
and so those limits were raised to 2.5xrms. In both cases, the Naturally-weighted cube with 
the standard CLEANing depth of 2.5xrms produced good results. 

For galaxies with all new data, we produced cubes at the observed spectral resolution and 
a second set with Hanning-smoothing applied. The Hanning-smoothed version is presented 
here. For cases where we were combining new data (not Hanning-smoothed) with archival 
data that were Hanning-smoothed, we shifted the new data in frequency space with the 
Hanning-smoothing option turned on (cvel). 

To remove the portion of each channel map that was only noise, we used the Natural- 
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weighted map cube to determine emission from the galaxy. We smoothed the original map 
cube to a beam of 25"x25" (convl) and blanked all pixels below 2-2.5rms of the smoothed 
cube (blank). We then made a second pass, blanking by hand any additional emission that 
was clearly noise (for example, features that do not appear in several consecutive spectral 
channels). This blanked cube served as the master cube that was used as a conditional 
transfer for the cubes with smaller beam size. We made flux-weighted moment maps from the 
resulting data cubes (xmom): integrated Hi, intensity-weighted velocity field, and intensity- 
weighted velocity dispersion. We corrected the integrated Hi map for attenuation by the 
primary beam (pbcor) and replaced blanks with values of zero (remag). 

A few galaxies were affected by Galactic emission. For VII Zw 403 Galactic emission is 
apparent in the C and D array data. The channels containing Galactic emission but no VII 
Zw 403 emission were blanked. However, Galactic emission in the channels that contain VII 
Zw 403 emission was not removed. For NGC 1569 Galactic emission is a serious contaminant 
over a large velocity range. See Johnson et al. (2012) for a discussion of how that galaxy 
was handled. 

We had particular problems imaging DDO 75. There were 8 uv data sets including AB 
(when the configuration is in transition from A to B), B, C, and D array configurations. Six 
of the data sets were archival and two were obtained under project AO 215. The moment 
maps have low level diagonal striping, and the integrated H I maps have individual blanked 
pixels ( "pinholes" ) where there should be flux. The increased noise from the striping caused 
the algorithm that produced the moment maps (xmom) to calculate unreasonable velocity 
values, and so those pixels were blanked in the moment map. We tried to identify the cause 
of the striping, but it appears at such a low level in the uv-data that we could not isolate and 
remove it. We concluded by cleaning to 3xrms since cleaning to 2xrms looked overcleaned 
and produced many more pinholes. 

Information about the ROBUST-weighted and Natural-weighted map cubes are given 
in Table [3} along with the channel separation of the observations and the final number 
of frequency channels in the cube. For each cube we give the beam full-width at half of 
maximum (FWHM) along the major and minor axes, the beam position angle, and the rms 
in the maps. The rms in a single channel varies from 0.46 to 1.1 mJy beam~^ in our sample. 
The pixel size of the cubes is 1.5", except for DDO 216 and SagDIG, which were observed 
only in C and D array configurations and were mapped with 3.5" pixels because of their 
significantly larger beam. The map sizes are usually 1024x1024 pixels. A few particularly 
large galaxies on the sky or galaxies with nearby Galactic emission were imaged at 2048x2048 
pixels. 
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4. Basic Analysis 

For each galaxy we have measured the total mass in H I, determined the H I flux- velocity 
proflle, and measured the Hi surface density proflle. These basic analysis products are 
presented here. To make these measurements, we used the Groningen Image Processing 
System (GIPSY, Vogelaar & Terlouw 2001). 

The cube channel map units are in Jy beam~^. This is related to brightness temperature 
in degrees Kelvin through 

5(Jy beam"^) = 1.65 x lO'^AaA^ TeiK) 

(Brinks, private communication), where Aa and A6 are the major and minor axis FWHM in 
arcseconds given in Table [3j The integrated H I map is obtained after converting the channel 
maps to units of K km s~^ using the above relation and summing along the velocity axis. 
The column density becomes 

Nniiatoms cm"^) = 1.823 x 10^^ ^ T(K) Al^(km s"^), 

which becomes 

NHiia.toms cm'^) = 1.105 x 10^^ S{Jj beam"^ m s"^)/(AaA(5) 

(Hibbard et al. 2001). The Hi mass proflle is obtained by integrating the flux in each channel 
map over the area containing Hi emission and dividing by the number of pixels per beam 
area to arrive at S expressed in Jy. The conversion to total mass is then 

Mhi{Mq) = 235.6 D(Mpc)2 ^ 5(Jy)Al^(km s^^). 

The factors for converting fluxes to masses for our ROBUST-weighted integrated H I maps are 
given in Table |4j 

4.1. HI masses 

To measure the total Hi mass, we summed the flux over the primary-beam corrected 
integrated Hi map made from both the Natural-weighted image cube and the ROBUST- 
weighted cube (flux). Since the integrated H I maps had been constructed from the blanked 
cubes and pixels without emission had values of zero, the flux was summed over the entire 
image. The flux was then converted to mass. The masses from the two maps agree well. 
The average of the absolute value of the differences between the masses measured from the 
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ROBUST-weighted maps and those measured from the Natural-weighted maps is 11%, with 
the ROBUST-weighted masses often higher than the Natural-weighted masses. 

The fluxes and masses from the ROBUST-weighted data are given in Table |4| and com- 
pared to single-dish measurements in Figure |3j The single-dish values are taken from Bot- 
tinelli et al. (1982), Schneider et al. (1990; DDO 101), Huchtmeier et al. (2003; DDO 53, 
DDO 216, F564-V3), and Springob et al. (2005; CVnIdwA). Our VIA Hi masses agree 
within 11% with those measured with single dish telescopes. However, there is some bias to- 
wards VLA interferometer masses being slightly lower than those measured with single-dish 
telescopes. One can see this in Figure |3] by comparing the best fit line to the data in red 
with the line of equal values in black. 

4.2. Velocity profiles and Single-dish comparison 

In general, we find good agreement between the VLA and single-dish measurements in 
the literature. We used spectra from the NASA Extragalactic Database, when available, 
and the telescopes that the single-dish measurements come from are given in Table [5) When 
comparing velocity profiles, we use our VLA robust maps as they generally recover more 
flux than the Natural-weighted maps, which were cleaned to 2.5xrms rather than 2xrms. 
The profiles are shown in Figure |4] and details of the comparison of individual galaxies are 
listed in Table [H 

Galaxies that show no significant difference between the interferometer and single- 
aperture measurements include CVnl dwA, DDO 43, DDO 46, DDO 47, DDO 52, DDO 
70, DDO 75, DDO 133, DDO 155, DDO 210, DDO 216, SagDIG, and VII Zw 403. Those 
that are very close in shape, but have 10% to 25% less emission in the VLA profiles include 
DDO 53, DDO 63, DDO 87, DDO 167, and DDO 187. Most of these differences fall within 
the range of measurement uncertainty, and the choices of baseline fits can also account for 
some of the discrepancies. Some galaxies (LGS3, M81dwA, NGC 3738, UGC 8508, Haro 36, 
Mrk 178) have rather noisy profiles, but generally agree in shape. 

Galactic Hi emission affects several of the single-dish profiles. For DDO 53, DDO 69, 
Sag DIG, and VII Zw 403, the addition is obvious and far enough away in velocity from the 
peak of the extragalactic emission that it introduces no confusion in the comparison. NGC 
1569, however, is heavily affected by Galactic H I, so that the single-dish and interferometric 
profiles are hard to compare. 

IC 10 and IC 1613 exhibit much more emission in the single-dish observations, with the 
VLA only recovering 40% to 50% of the peak single-dish flux. The VLA could be missing 
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a significant amount of emission simply because the galaxies have large angular extents. 
Interferometers are inherently less sensitive to emission filling a substantial fraction of the 
primary beam. DDO 50, DDO 69, DDO 126, DDO 165, NGC 4163, and NGC 4214 are 
other galaxies that show more modest examples of this behavior. The high-velocity gas in 
Haro 29 matches quite well with single-dish observations, but at the low-velocity side of the 
profile, the single-dish observations have a broad wing of low brightness material that does 
not show up in the VLA profile. 

For galaxies that occupy a larger area on the sky than the half-power beam width of a 
single-dish instrument, we expect the VLA, with its relatively large HPBW, to recover more 
extended fiux. This can show up in the form of enhanced "shoulders" or asymmetric excesses 
in the VLA data. Galaxies in this situation are DDO 101, DDO 154, DDO 168, F564-V3, 
NGC 2366, and WLM. Two single-dish profiles are presented for DDO 154 because there is 
some disagreement between the published single-dish data. 

4.3. Surface mass densities 

To construct the H I mass surface density profiles, we measured the fiux in annuli on the 
velocity- integrated Hi maps (ellint). We did this using the optically-defined parameters 
of center position, position angle of the major axis, and inclination of the galaxy. These 
parameters were determined from the optical y-band morphology (Hunter & Elmegreen 
2006), and are given in Table |6} Hi kinematic parameters are in the process of being 
measured as the velocity fields are determined (Oh et al., in preparation), and eventually 
surface density profiles will be made with the kinematic parameters. The optical inclination 
assumes an intrinsic minor-to-major axis ratio h/a of 0.3 (Hodge & Hitchcock 1966, van den 
Bergh 1988). The width of the annuli and step size for the radius of each annulus for each 
galaxy was chosen to be approximately that of the FWHM of the beam, and the profiles 
were integrated from the center of the galaxy outward until we ran out of image. The surface 
density profiles for the ROBUST-weighted integrated H I maps are shown in Figure [5j 

5. Public data products 

The Hi and ancillary data are being made available to the public through NRAO. The 
Hi data include the Naturally-weighted and ROBUST-weighted cubes and moment maps 
(velocity-integrated H I maps H I, intensity-weighted velocity field, intensity-weighted veloc- 
ity dispersion). The ancillary data include FUV, NUV, UBVJHK, and Ha images, as 
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available for each galaxy. In addition to the images and cubes, we will include basic analysis 
products in tabular and graphical form as they become available. These currently include 
the velocity-integrated flux profiles and H I mass surface density profiles. 



The NRAO web site is http://science.nrao.edu/science/surveys/littlethings, and the 
data can be obtained there now. The ancillary data are also available at the Lowell Obser- 



vatory web site http://www.lowell.edu/users/dah/littlethings/index.html We believe these 
data will be useful to the community for a wide variety of investigations. 
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A. Calibration and Combination of Data 

The calibration steps followed broadly standard procedures; here we highlight those 
steps that needed particular attention and our "recipes" are available online. In what follows, 
CHO refers to an average of the inner 75% of the frequency channels observed; LINE refers 
to the spectral line data set. 

Data were loaded into AIPS using fillm with the option to store channel/IF dependent 
weights with the visibilities, creating a CHO and LINE file. The CHO data were inspected 
for interference and other obvious problems using tvflg. The visibilities were updated with 
baseline corrections from vlant. The EVLA-EVLA baselines were flagged in their entirety. 



- 21 - 



We removed 8% of the channels, 4% on either side of the band, from the LINE data 
using UVCOP because the noise in the first channel would percolate through to the bandpass 
calibration. We subsequently proceeded with the bandpass calibration. Contrary to what 
was customary when reducing VLA data, we did the bandpass correction as one of the first 
steps, rather than as one of the last. This was to remove VLA-EVLA closure errors due to 
non-matched bandpass shapes between VLA and EVLA antennas. 

We recreated a new CHO using AVSPC on the LINE data and applying the bandpass 
correction. It is this CHO that we then further calibrate in a fairly standard manner. We ran 
a number of diagnostic tests to check the quality of the calibration and to detect possible 
problems, notably those related to the VLA-EVLA transition period Quality control checks 
included using POSSM to examine the band-pass calibration done prior to constructing a new 
Channel (CHO) (to examine the band-pass calibration table, to apply the calibration to 
the secondary calibrator and examine individual baselines, and to apply to the secondary 
calibrator as a vector average of all data), snplt to inspect the amplitude and phase cal- 
ibration (SN) table, listr to look at the SN table itself and check for phase or amplitude 
jumps, ANBPL to check the antenna-based data weights, tvflg to examine the uv-data as 
a function of baseline and time, uvplt and wiper to plot the amplitude and phase against 
uv-distance for the calibrators and galaxy with the CHO and with the line data, IMAGR to 
map the calibrators and check their fluxes and map the source with the CHO and line data, 
POSSM to check the vector average of all of the data for each calibrator, and computing the 
actual and expected S/N. 

With few exceptions (see Table [s]) the LITTLE THINGS galaxies were observed in 
three different array configurations, B, C, and D, that needed to be combined into one data 
set. Furthermore, multiple observing runs contributed to the data taken in each of the 
configurations. Some of the data came from the NRAO archive, whereas the rest were newly 
observed. The archival data were Hanning-smoothed whereas the new data were not. 

Once calibrated, we ran task CVEL to correct the frequencies from scan to scan as well 
as to align all observations of each galaxy in velocity. If newly observed data were to be 
combined with archival observations, Hanning-smoothing was applied to the data as well, 
again as part of task CVEL. Although this task applies Hanning-smoothing, it does not 
remove every other channel, which is what used to be the case with the VLA correlator 
when specifying the Hanning-smoothing option. To accomplish this we used task UVDEC. 
Once the data sets had the same frequency resolution and sampling, they were aligned in 
velocity and, where required, precessed to the same J2000 epoch using UVFIX. We used task 
DBCON to merge them into the final dataset. 

We used UVLSF, and a linear baseline, to extract the continuum based on the location 
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of channels free from line emission as identified in a dirty image cube of all the combined 
data. 



B. Imaging 

Interferometers measure the Fourier transform of the sky intensity distribution. How- 
ever, interferometric data are affected by incomplete coverage of the aperture plane: missing 
zero-spacing flux and information on short-baselines, missing or deleted baselines, missing 
or deleted hour angle ranges. To remove the effects this incomplete coverage has on the 
resulting images, deconvolution algorithms were developed, notably the method known as 
CLEAN (Hogbom 1974). Although powerful, CLEAN has its limitations, some of which we 
mention briefly here before we describe the approach we took. 

Standard CLEAN decomposes an image into clean components (basically Dirac delta 
functions at grid point locations) and a residual map. The clean components convolved with 
the original (also known as dirty) beam, a process known as restoring, and put back on to 
the residual map would give the original (dirty) image back. If instead the restoration is 
done with a Gaussian beam that matches the resolution of the dirty beam, a cleaned map 
results. Usually the flux in the cleaned map is determined under the assumption that the 
clean beam applies, which is only true if all flux registered by the interferometer was cleaned. 
Because noise in the map makes it often difficult to clean down to the noise level, the residual 
map can contain appreciable residual signal, especially when cleaning extended structures 
(Jorsater & van Moorsel 1995). The flux is calculated correctly only for the clean components 
retrieved from the map, for the residual map it is overestimated by a factor equal to the 
ratio between the dirty and clean beam (for more details, see Walter et al. 2008). For this 
reason, a standard clean using AIPS task IMAGR can overestimate the flux by as much as 
50% as shown in Fig. |6] (bottom left panel, red line). Rescaling the flux in the residual map 
can be used to correct for this (Fig. |6| bottom left panel, green line), but does not deal with 
the fact that a substantial amount of flux has not been retrieved, and hence cleaned. This 
manifests itself by the emission sitting within a depression in the map (commonly referred 
to as the negative bowl), and in Fig. [6] (bottom left panel), by the subsequent drop after 
reaching a maximum when integrating a map out to ever larger radii. 

Extended, low-level emission is a major contributor to the integrated flux of extended 
sources, even if the brightness in each pixel is at low signal-to-noise, which means that 
CLEAN would have to go to very low levels to recover all signal. That is difficult as CLEAN 
tends to diverge when reaching levels below typically 2xrms. To address this, Wakker & 
Schwarz (1988) developed an algorithm that solves sequentially, in the image plane, for 



- 23 - 



structures at increasing scale. Their algorithm formalized the concept introduced by Brinks 
& Shane (1984) under the name of Multi-resolution CLEAN, who ran CLEAN on an image 
at full resolution down to a set limit, after which the residual image and dirty beam were 
convolved to lower resolution and cleaned further, and so on, cleaning at 3 or 4 resolutions 
in total. 

An extension of the Wakker & Schwarz (1988) algorithm was further developed by 
Cornwcll as Multi-Scale CLEAN (Cornwell et al. 1999; Cornwell 2008) and implemented 
in AlPS by Grcisen as a Multi-resolution CLEAN within AlPS task IMAGR (Creisen et 
al. 2009). An early implementation within the Common Astronomy Software Applications 
(CASA) package was tested on galaxies that are part of THINGS (Rich et al. 2008). We will 
refer to both of these implementations as Multi-Scale CLEAN (or MS-CLEAN for short). 

MS-CLEAN deals with extended structure by using a priori knowledge of spatial scales 
expected to be present in the image. When smoothing, the ratio of the extended source 
structure to beam size improves, which means that what was extended in the original reso- 
lution becomes small scale in a smoothed map. Moreover, the signal to noise in a smoothed 
map improves as long as the structure in the map remains resolved, allowing a much deeper 
clean and therefore ensuring that most flux is recovered. For information on all scales to 
be retrieved in the AlPS implementation of MS-CLEAN, the user not only can choose the 
widths of circular Gaussians that will be used to make beam and data images at each scale, 
but also the relative importance given to small versus large scales. 

The AlPS implementation of MS-CLEAN uses the multi-facets approach available in 
IMAGR for wide-field imaging. Whereas for wide-field imaging iMAGR creates a map made 
of adjacent facets, in MS-CLEAN mode IMAGR creates a map out of overlapping facets over 
the same area of the sky, each facet representing a different resolution scale instead. At each 
major cleaning cycle, it selects one of the angular scales, re-images that facet with the current 
residual — data, finds clean components for that facet using minor cycles and continues 
until it reaches one of the stopping criteria, prompting clean to go through the next major 
cycle. The selection of which resolution should be cleaned is done by estabhshing which 
resolution has the highest peak flux. The algorithm stops when the flux level in all facets 
(at each of the resolution scales) is below the user-defined flux cut-offs or the number of 
components cleaned has reached a set limit. Then, the clean components for all resolutions, 
each with their corresponding clean beam, are restored to the image made at the highest 
resolution (Greisen et al. 2009). 

It should be stressed that the parameter choice should be carefully tuned to the data 
set to be imaged. The algorithm available under the IMAGR task in AlPS requires a number 
of input parameters. The first two of them are ngaus and WGAUS representing the number 
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of scales to clean simultaneously and the approximate width of the convolution kernel used, 
in arcseconds. We applied the algorithm to images with data collected in B, C, and D- 
configurations, i.e., with a substantial range of angular scales (spatial frequencies). After 
several tests (see Fig. [6| top left panel) we decided to use four kernels with widths of 0" (i.e., 
native resolution), 15", 45", and 135". 

FGAUS allows the user to impose a stopping criterion for each resolution. We established 
our rms noise threshold by measuring the noise at each of the four resolutions of a dirty, line- 
free channel. Consequently, we have four different noise levels, one for each field. The four 
cut-off levels to be set in FGAUS are not independent, being different smoothed versions 
of the same dataset. We investigated different flux cut-offs, trying both deeper cleaning, 
down to Ixrms and 1.5 x rms and more superficial cleaning to 3xrms and 5xrms. We plot 
the integrated flux as a function of radius in Fig. [6] (top right panel). We find that cleaning 
down to 3 or 5xrms is not sufficient as we are not retrieving all the flux and evidence for 
a "negative bowl" is seen at large radii. Cleaning down to 1 or 1.5xrms does not retrieve 
more flux than cleaning down to 2xrms; however, clean artifacts start to become noticeable 
in the final image. This implies we do not need to go any deeper than 2xrms to retrieve all 
the observed flux. We therefore decided for LITTLE THINGS on cut-off levels of 2xrms for 
ROBUST-weighted data cubes and 2.5a for the Natural-weighted cubes. 

Naively one would always look for the peak residual at any resolution, and choose the 
next clean component based on the highest value found. However, since the lower resolutions 
dominate in terms of flux, the algorithm would be biased towards the largest scales (Greisen 
et al. 2009). To balance this, peak fluxes can be weighted by a factor of 1/ [(field beam 
area) /(minimum beam area)]". Here a provides control over which resolution is going to 
be chosen next. If a = 0, the peak fluxes at each resolution are unaltered; if a = 1, the 
peak flux at the highest resolution will be the true peak flux, while the peak flux at all other 
resolutions will be down weighted. In the AIPS implementation of MS-CLEAN, a is set by 
the parameter imagrprm(II). 

To illustrate how this biasing works we plot the choice of resolution in each major cycle 
when cleaning one channel (Fig. [?]). We ran tests on the choice of a and we found that for 
all a, the cleaning process converges. Also, the total retrieved flux is very similar. We found 
that a value of a = 0.2 strikes the right compromise, with the clean alternating among the 
different spatial scales and thus finding emission across a range of angular scales rather than 
being fixated on one particular value. The effect of different values for a on the integrated 
flux is shown in Fig. |6] for both robust and Natural-weighted data. 

The differences between a = 0.2 or 0.4 become more noticeable when looking at the final 
and residual maps resulting from applying MS-CLEAN to a line channel for two different 



- 25 - 



sources, DDO 168 and DDO 133 (see Figs. [8]and|9]). It is clear that a particular choice for 
a depends on the complexity of the source to be imaged. Although there is little difference 
in the final maps for DDO 168, in the case of DDO 133 there are signs of over-cleaning for 
a = 0.4. 

As an aside, it is interesting to note to what extent the residual scaled standard CLEAN 
follows the Natural-weighted MS-CLEAN results, even at larger radii, which was not the case 
for the ROBUST=0 MS-CLEAN case (see Fig. |6| bottom right panel). This is not an effect of 
MS-CLEAN, but rather a consequence of the different density of visibilities near the center of 
the Mf— plane and the different manners this is dealt with in the weighting scheme employed. 
Natural versus robust. 

Other than those described above, all further parameters were left at their default 
values. If we compare MS-CLEAN and traditional CLEANing methods we find that in our 
test the standard CLEAN has a peak flux which is 40% higher than the MS-CLEAN or 
residual scaled CLEAN (see Fig. [6| bottom left panel). This shows that MS- CLEAN, unlike 
standard CLEAN, properly addresses the problem related to determining the total flux. 
At the same time MS-CLEAN recovers all measured flux much better than the standard 
CLEAN, reducing considerably the negative bowl problem in extended sources. As a caveat, 
MS-CLEAN does not recover flux that was not measured at spacings shorter than the shortest 
spacing partaking in the observation. Therefore, any emission on scales more extended than 
those present in the observation will be missing. 



we 



As a final check, we looked into the noise characteristics of MS-CLEAN. In Fig. 10 
present noise histograms for different choices of a corresponding to the cleaned maps at full 
resolution, and the 15" and 45" residual maps and compare these with noise histograms on 
a line free channel at the same scales. We see that MS-CLEAN does not affect the nature 
of the noise in any of the maps shown (the tail to positive intensities in rows 2, 3, and 4, 
reflects genuine emission). 



C. The Data 



are 



We present the basic data for each galaxy in Figures 11 to 90. Figures 11 and 12 
shown here, and the rest are available on-line. Each galaxy has two pages: the channel maps 
and a combination of images. The requirement for the channel maps was that each galaxy 
fit on one page, and so some have binned pixels and none show every channel. However, we 
show the full range of channels with emission. On the image page, the false color image in the 
upper left combines the integrated H I map (red) , V (green) , and FUV (blue) and allows one 
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to see the relationship of the gas and stars. The other images are the first three moment maps 
from the Naturally-weighted cube: velocity-integrated H I, fiux-weighted velocity field, and 
flux-weighted velocity dispersion for each galaxy. All of the fiux-weighted velocity dispersion 
maps (moment two) were made with a scale of (dark blue) to 15 km s~^ (red). 
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Fig. 1. — Properties of the LITTLE THINGS sample (new VLA observations [hashed +45°] 
and archive data [hashed —45°]) compared to the entire optical survey of Hunter & Elmegreen 
(2006; no hash). The LITTLE THINGS sample covers the range of parameters of the full 
survey. SFR^i is the star formation rate, based on Ha emission, normalized to TcRj^, where 
Rd is the disk scale length determined from a V-band image. A star formation rate per 
unit area of is plotted as a log of —6. The parameter fi^ is the central surface brightness 
determined from a fit to the y-band surface brightness profile. The distance dependent 
parameters My and Rd use the distances given in Table [T]for the LITTLE THINGS sample 
and the distances given by Hunter & Elmegreen (2006) for the rest of the galaxies in the 
larger optical survey sample. 
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Fig. 2. — Noise level as a function of channel number for data obtained in B-configuration 
on DDO 43. The total bandwidth was 1.56 MHz split into 256 frequency channels. The 
observations were made with 26 antennas, 12 of which were equipped with the new EVLA 
front-ends, during the VLA-EVLA transition period. We show here different combinations 
of baselines: combination of 66 EVLA-EVLA baselines (navy crosses); 91 VLA-VLA base- 
lines (green triangles); the sum of 168 VLA-EVLA baselines (open magenta squares); the 
combination of the VLA-VLA plus VLA-EVLA baselines (red stars); and finally all 325 
baselines (open blue diamonds). Continuum emission was aliased, with phase information 
scrambled making it resemble noise, the effect being strongest in channels closest to the 
baseband edge. This is most noticeable on EVLA-EVLA baselines, and to a lesser extent on 
EVLA-VLA baselines; VLA-VLA baselines were not affected. 
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Fig. 3. — Comparison of our integrated H I masses with those from single-dish radio telescope 
observations. Both measurements are relative to the same distance. Our Hi masses are 
measured on the VLA integrated Hi ROBUST-weighted map. The solid black line marks 
equal masses, and the red dashed line is the best fit to the data. The average difference 
between the masses measured from the ROBUST-weighted maps and those measured from 
the Natural- weighted maps is 11%. The best fit line lies slightly below the line of equal 
masses, indicating that the VLA has systematically missed some flux compared to single- 
dish values. 
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Fig. 4 

of the 
Single- 



:. — Intensity- velocity plots for each galaxy made from summing the flux in each channel 
i Hanning-smoothed Naturally-weighted [red) and ROBUST- weighted [blue) Hi cubes. 
!-dish profiles are also shown {black; see Table [s] for references). 
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Figure 4 (continued) 
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Fig. 5. — Azimuthally-averaged Hi surface density as a function of distance from the center 
of the galaxy. These are measured from the ROBUST-weighted velocity-integrated H I maps, 
using parameters — center position, position angle of the major axis, and inclination of the 
galaxy — determined from the optical \^-band morphology (see Table [6] Hunter & Elmegreen 
2006). The width of the annuli and step size for the radius of each annulus was chosen to 
be approximately that of the FWHM of the VLA beam. 
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;ure 5 (continued) 
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Fig. 6. — Integrated Flux density versus radius in DDO 168, channel 60, B+C+D- 
configurations combined. Top Left: Results obtained using MS-CLEAN with a = 0.2, robust 
= 0, and testing different combinations of angular scales (red: 0, 15", 45", 135"; blue: 0, 15", 
45", 135", 407"; cyan: 0, 20", 60", 180"; magenta: 0, 40", 80", 120"). Top Right: Results 
obtained using MS-CLEAN with a = 0.2, robust = 0, and testing variable flux cutoffs (red: 
la; blue: 1.5o"; cyan: 2o"; magenta: 3cr; navy: 5a). Bottom Left: Results obtained by clean- 
ing down to 2a, robust = 0, and testing dependence on a (red: standard CLEAN, green: 
residual scaled CLEAN, blue: MS-CLEAN and a = 0, cyan: a = 0.2, magenta: a = 0.4, 
orange: a = 0.5, navy: a = 0.7, brown: a = 1). Bottom Right: As Bottom Left, but with 
natural weight. 
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Fig. 7. — Testing the a parameter on DDO 168, channel 60, B+C+D-configurations 
combined. Wc plot the choice of resolution scale in each major cycle for a — 
0,0.2,0.4,0.6,0.8,1.0. The parameter a provides control over which resolution is chosen 
next. If q; = 0, the peak fluxes at each resolution are unaltered; if « = 1, the peak flux at 
the highest resolution will be the true peak flux, while the peak flux at all other resolutions 
will be down weighted. We found that a value of a = 0.2 strikes the right compromise, with 
the clean alternating among the different spatial scales and thus finding emission across a 
range of angular scales rather than being fixated on one particular value. 
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Fig. 8. — DDO 168, channel 60, B+C+D-configurations combined: Comparison between 
2xrms cut-off level, robust 0.5, clean maps (top) and residual maps for 15" {second row), 
45" {third row) and 135" {fourth row) for a — 0.2 {middle column) and a — 0.4 {right 
column) . The beam size is shown in the bottom left corner of each figure of the middle and 
right columns. In the first column we show as reference the uncleaned map for each field. 
All panels of the figure use the same grey scale, shown on the right hand side. 
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Fig. 9. — DDO 133 channel 56, B+C+D-configurations combined: Comparison between 
2xrms cut-off level, robust 0.5, clean maps (top) and residual maps for 15" {second row), 
45" {third row) and 135" {fourth row) for a = 0.2 {middle column) and a = 0.4 {right 
column). The beam size is shown in the bottom left corner of each figure of the middle and 
right columns. In the left column we show as reference the uncleaned map for each field. All 
panels of the figure use the same grey scale, shown on the right hand side. 
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Fig. 10. — Histogram of the intensities in a single channel of DD0168, B+C+D- 
configurations. The row at the top is for reference and represents the distribution in a 
line-free channel, channel 87, robust — 0.5, and without any cleaning applied. The first 
column corresponds to the full resolution maps, second column to the 15" residual map and 
the third column to the 45" residual map. Rows 2, 3, and 4 are histograms of a line chan- 
nel, channel 60, cleaned using MS-CLEAN with different values for a, from top to bottom: 
a = 0, a = 0.2, and a = 0.4. Note that the first column corresponds to the final (cleaned 
and restored) map which explains the tail to high intensities; the remaining results are based 
on residual maps. 
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Fig. 11. — Channel maps for CVnIdwA made from the Hanning-smoothed Naturally- 
weighted Hl-hne cube. The synthesized beam is given in the left bottom corner of the 
top left panel. In each panel we give the observed radial velocity in km s~^. 



- 45 - 




Fig. 12. — Basic data for CVnIdwA. Hi data are made from the Hanning-smoothed 
Naturally-weighted cube. Upper left: False color picture combining Hi (red), V^-band (green), 
and FUV (blue). Upper right: Integrated Hi. The color wedge is in units of Jy beam~^ m 
s~^. Bottom left: Intensity- weighted velocity field. Bottom right: Intensity-weighted velocity 
dispersion. The color wedge for both bottom panels is in km s~^. 











cci 

























+ 












t) a 

cd - 
1 

M ^ 

o 


















0) 1 

« s 

cd 

^ 7 
^ I. 
























2- 






a 












E 










'Ji 

'x 















u 

Q ^ 












£ 
1=1 



















^. 














46 



q 


q 










i-H i-H O (T) CT) ^ CO 

i-H i-H i-H O O O O 




S S i£ 5£ 






lO 


LO 

^ o 






q 


q 


q 




q 






d 


d 


d 


d 


d 




O O O O O O O 


X 






C: O 




o 


q~ 


CO 


d 


d 


d 


d 


d 


d 


d d 


-H 


-H 


-H 


-H 


-H 




^ ^ ^ ^ ^ ^ ^ 




^ ^ ^ ^ 




^ -H 


-H 


-H -H 






-H 


-H 


-H 


-H 


-H 


-H 




CO 


CO 




oo 


t- 




^3 ^ ^ !!?■!!?■ OO t- 




1.0 t> CD 








CM CD 


















CO a 




00 


00 








t~ t- t- t- t- t- 00 




00 




00 1> 




00 1> 






00 




00 




00 


00 





a>_ c:>_ c:>_ c:>_ c:>_ c:>_ a a a a a a a <o_ <o_ <o_ 
Sdididididi<6<6<6<6<6<6<6<6<6<6 
-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H 



o_ <o_ <o_ a a a 
<6 <6 6) (6 (6 (6 
-H -H -H -H -H -H 



oooooooooo 



oooooooooo 
-H-H-H-H-H-H-H-H-H-H 

OJOOtOiHtOM'^MiN'-i 



I I I I I I I I I I I 



I I I I I 



qqqqqqqqqqqqqqqqqqqqq 
ddddddddddddddddddddd d 
-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H-H :-H 



-H -H 



Q Q Oi^ Oi^ a <o^ <o^ 
'■ 6> <6 <6 <6 

i-H-H-H-H-H^+l+l 



OOOiOOOOOtN' 



ooooooooooooooooooooooooooooooooooooo 



OqCOCOCOlOCOtO'H'H'H'HCSCOCOC^COi-lXLOiH'H.H.HCO'HCNCNOC^l^iHCO'HCOCO'HCO 

t-hOOOOt-iOOOOOiHOCCDOOOOOOOOOOOOOOOOOOOOOO 

ddddddddddddddddddddddddddddddddddddd 
-H-H-H-H-H-H^^^^^^^^^^^^^^-H-H-H-H-H-H-H-H-H-H-H-H-H-H^^^ 

t:^I^^I^OOCSOOO^OOC^'^^^~'^C^ln!DeO(^lOO^-'*COOOOCO^DOO^DOOI^>OI^CO^-t^ 

lO^IHcO'Hco^~Cl'H"sl^c^coolooMl^)'Hc^cooO'HIHlOl0^lOcs^^cccc^-^^^-lO^ 
odiHiHiHiHodddd^dd'Hdddddddddddddd^ 



iNC0iH,-i<N<NMM^^'-i(N'-iO(NO(N'H,H-^ 



oi d 



tN O t- 

^ ^ CD 

IN iH 



'#T-lt~lOCD'*OOOOt- 
(NiO^iodiOCO^ 



m q q cr; 
CO in lo 



CN in CO in CN 



ILN iiJ ixj lij lij T-H -qi iij iij -rj" -rj" -rj" UJ C'J UJ (LN CJ t^J ^ l-^ W i^N 

tHiHiHiHiHiHi-Ii-Ii-Ii-Ii-Ii-Ii-I.H.H.H.H.H.H.HiHiHiHiHiHiHX'iHiHiHiHiHiHiHi-1 



(NOOi-i^^DtMintD^ 
oddt^^t^lOlNCO^ 



MiH'^i-li-llOCOl^ 



.^GOGOiH^^-^^^^^^^COiHiH^^. 



tDXIHl^|■^^•!^ooioococo^-■^cKlO^-c^coc^coc^Ol'H^-^-^-^-lo■^■^OlOlOl 
co^-^tdlOco3cocod'-^'-^^-^cdTl'cococ^■^■^■^c^d'^ 



ni 

00 



O d 



O o ■ 

O O . 

D O , 

in CD 

00 o 

in CD 



(M (M (M 

o o o 



Hi Hi Hi 



in CD 

O 

CO CO 

^ (M 

O O 

O O 

D Oh 



in. 

^D 

m O " 
in in 

S S oo' 

^ ^ CD 
O 00 
lO G 00 

o "3 

Sao 

O o 
O 

o 



00 ^ 



:o "D 
o 

CO 

< 



ro 

O ro . 



O 2 

oSro-§ 
O J ro CO 

^ . ^ 00 

o o -"so 

O O S ^ N o 
p p o O H p 



ro 00 ( 
in ro 
in CO ( 



00 



O 00 > ^ 

N N 5 _ 

CO -H Q fh 

o o .2 Q ° 

^ ^ CO 

^ ^ ^ n, in 



O ■ 
N D 
2 D 

O " 

o S 



o o 
o o 



o 

O o 

ro o 

CO ^ 

00 P 

CO . 

^ ^ 
C31 

U 

o 
o 



s; S a a ri <; 



o o 
•p o 



CO CO 

o ro 

CO o 

^ in 

O O 

o o 



S o 

CO CQ 



CO o 

o o 

^ 00 
CO 



o 



lij w 
a. „ o 



■ o 

rQ c 

^ ^ m 

m CD 

' O) in 

■ CO CD 

' O O 

; o o 



' o m 

' ^ °£ 

^ in 

^ CN CO 

; o o 

! a a 



O m 

t-' O 

CD O 

5 D 

O in" 

g I 

S o 

00 o 

CO CX 
O 



^ CO i: 

ffl O CD O 

ro" CO ^ ^ 

c^ -^^ - 

C^ CC M ^ 

O o in c-i 

Q ^ 3 o 

- Q 



g S u 

O) CD O 

O cu 

O 
(X 



CD 





CD 




O M 


CO CO ro 


^ CO CO ^ 1C5 

oinb-ocicoinin 


1C5 
CO 


CO 


« t~ o 

CO 00 tH 


^ 






in in 


in CD CD 


b-b-OO^^^^^ 






^ ^ 

















00000000 









Q 


Q 


Q 


Q Q 


Q Q Q 


QQQQQQQQ 


Q 


Q 


Q Q Q 


Q 


Q 


Q 


Q Q 


Q Q Q 


QQQQQQQQ 


Q 


Q 


Q Q Q 



Ol !0 W P3 ^ (N OO 

jcocomtOi-itN o 

' "^lOCOt-i-llMOOrKin 
T-\ ^ & i-l(MCO'#'#CD^0O 

!riri 0^000000 ^Ol^ 



^ CN ^ 



O O O C 

^ ^ ^ ^ 

^ b- t> 

t-^ o6 t-^ t> 



9 9 9 9 

O: O: O: d:i 

b- m CD t- 

o m CD CD 



9 9 9 9 

<6 <6 <6 <6 

+i +1 +1 -H 



o o o c 



o o o o 

<6 <6 <6 <6 

^ ^ ^ 4^ 

CO CM 

CN CD CO in 

S S S d:^ 



CD i-H CO 

^ m ^ ^ 



oo CO C7i ^ 

m d CO ^ 



55 in ^ m 

^ t- CD CO 

5 O O O 
^000 

D D D 

°o . . . 

^ ^ ^ CO 

^ CN oo 00 

I— I CD CN 

n CO m Lo 

Q ^ CO CO 

^0 

g H H H 

CO 



O O _^ N 
ni ni iH H 

K K S > 



cd 3 



C5 



Pi 



01 ^ 

J3 Pi 

s o 

a; " 



01 > 

IS 

cq 
a 



S -a 



s 



S Pi 5 



CN CD 

s ° 

o o 

CM IN 



Q £ 



o m 3 



^§1 



s o s 

— o 

n3 O 

+j (N . 

a> o 

CM -H ^ 



c^' 's 
c 

S o 

■H Ol _ 

3 " 

O 

S 

J ^ 

4J CO 

3 2 S 



e S 
S S 



« s V . 

4J ^ 00 

0) ^ ^ GO 

- Ol 

C CO ^ 

(1) ni ^ 

-id 

.£ « « cs 

O « " « 

" g S » 

O tN 



CN ISl 

O 0) 0) 

: W J r3 

CD <D W O 

5 K ^ ^ 



^ -4 CO 

' c6 00 



o W D H 



2 ro " S 

^ ni in 



-J CO ^ 



■5 N « 

2 S-' 



-48- 



3 



CO W CQ ^ 

:3 ^ 2 13 3 
H H H H K 



OJ CD 



to 
o 
o 



T3 a, 

CD r3 

o o g 

T3 T) 

a> oj g 

^ ^ 

cS <i3 

4J +J CO 

cS cS ^ 

Q Q O 



-o 
-% 

X 

bC CD 

O 



3 
-a 



3 



3 





c 













+J 


o 






i 




a 








a 









ati 












a 
















CO 




^ 






CD 






Sh 






S-i 


ec 




u 




a> 


q; 




o 


u 


o 


cS 










Tr 



CD q; 
o o 



^ 13 O 



H H m 



<A 

O 

6 



ce 
Q 



CO 

a ci. 

5J CO 

;=! > 

< < 
> 



O 

CD 



bO bO 
c^ cci 

a a 



M CO X! 

CD CD O 

a a ^ o o K 

I ^ " 



Q 



^ Pi a 
_o o <; 

O 



3 

Q 



-49- 



CM 



hH iT W hH 

S <: <: 

lO K >H lO 

S <: 
o 
<; 



►J o 
<; <; 



K K 

g <: <: 

<C o 01 

^ 

cn cn 

<; <; 



CM 



cncncncNCN'^cNi-ii-ii-ia;a;a;c3itototo,H,H,HOcMcM 



O n:^ « O ■ 

^ CO in rr r- , 

_ ^ lo t> c/3 *^ ; 

OClHClHClHffiffiffiffi^^^CUCUCU>Hjjco'^KH^'l 



~c^^oc^^c^^t^^I-^c^^c^^c^^ 

^ hd hd hd hC 



O < 



° -g ° 



<; O « : 



N " S S 2 ^ f- 



' 8 >2 -s a 

rP 2 6, << 



' < 

) ^ O) 01 r-i c-i 



i P S « £ o 



> Q & 



,CO(M1^c-i^01go"(^cO 

'^►^PSP.^^g^P 



c o o ^ as' 

CO o 



^ ^ CO ^ ^ — ■ - - 



IS S ^ 

s g g ^ 

- d" M go 

^ ^ ^ _£) 00 

; in ^ CD cri _ 



S h 2 

B S o 

3 ° « 

^ of O 

00 ^ 

- 00 

S o 



o 

" 5 !^ 

w ^ ^ 

1 « I 

-> s 

' S 2 00 



m^opmopmop^ogpnopnBOpnopnopmopmop^mopmopmBOPMo 









o 

lO 


lO 


CO 


CO 

to 


01 

CD 


o 
^- 


in 
t- 


l> 

00 


o 




DDO 



P 
P 




p 
p 



P 
P 



P 
P 




p 
p 




p 
p 


P 
P 


P 
P 


P 
P 



P 
P 



P 
P 




p 
p 



- 50 - 



ni X 
0) 



CO I> '^i q (N CO 
i-i i-l 2 to 00 to 



l> lO iH O CO 



rjt 00 ■ 



Oi <^ ( 



lOOOOiHtOtOiHiHlOOb-lNCO' 



I (O 1^ lO 1^ (N lO (N 



iH iH M IN CO 



^(0(Nt-'Hcocotoio(NcoioiNoicoco^;i|ioo^inco'*'*cotoinco 



(N(N(N(N^OOO,-i(Nt- 



(N (N 
01 01 

t>t>t>int>|i31il(N,Ha)(^i'-H0i'^'^i'^i^^'^^^^^ 

CNCNCN'HCN^^lf3l>.T^lOt^-^tN1010"^QOLOir:;000 

a;a;a;a;a;ajc::c::cr:cr:cr:cNc:^'\'\t>^cot-^coa)(:-t~cDcocoM^^^^, ^ ■ — .i- 

■ ' 



b- CT) Pi |i] CM 
^ ^ ^ 9 --^ 



S 2 « " 



00 S 



< 2 



O g O g 00 



O 3 ni ni c 



N e ^ 



O J3 = 



00 g 



" 1- 

< 00 N 



S ^ ° 3 



oo °0 
g O 



0)' ^ 
^ 00 o ' 



2 1h 



2 s <o-<: 



t- s 



lO O ^ to 



Z d oo < 



^ ■ 



fl 00 



g>3 S 

g M 



Q Q m 



o 

CO 


00 








of 

(N 




in 




bO 
=! 


Dec 




in 




Apr 


to 


u 

OJ 

Q 


o 




01 




a 

nJ 


> 



CO 

o 


u 
o 


(N 
O 


u 




to" 




Z 


OO' 


Q 


> 


Q 


^ ^ 




to lO 


00 

o 




CO 




tJ 


01 01 




Ol 


)h 


o 





o 


< 


0) 

Q 








(X 

< 




z 




01 


in 

Ol 








CO 

o 




o 





^pmoQMOQMsoQnBOQnsopmopmBciQmo 
Q o o o o 



QopMOQ^Mopnopnopmop 



-51 - 



o 



O iH lO lO 

;6 lo (N i6 d iH ! 



IN "vf 



a 00 

o L-^ o o lo t-- t-- w W 00 00 00 

OC^C^OOO(N(MCNiH(M 







00 


00 


lo ^- 


s 


X 


CO 


in CO 


a 




(M 


lO 
CM 


^ Ol 
Ol G> 


< 


< 


t- 


01 lO 
CO 










O X 


in 




X 




<; 


CO 






< < 


CM 


CM 


< 


< < 



o o 
< < 



o o 



O 



S « " S N 



CO fcj '-n 



§gS8§ 



<D 50 o > 



GO 00 

o o 



< 

H 00 

- o — 



01 » 01 



3 



o o 



IN < ^ 

§§§ 



^ o 

-V, °0 



a ■ 

<: ^ Q o>: , 

(N 52 I^- ^ ■ 

a <y> a V 



:S o 



01 



"5 CO 



CP O Q CQ 



OQOQCQ^OQCQOQ 



1569 


2366 


3738 


4163 


D 



NGC 


NGC 


NGC 



7-K lO 

Q O 



CI t- 01 



t- 1^ CO 1^ O 



lO CM 



CM O 



' QO 



CO 



00 CO CO 

" g ^ " 00 



- o < 

S ^ 00 00 



° o 5 o ^ o 



Ol" o" m <C 2; 
oo w 7" to t- 
^ 3 O 01 01 

00- tf S 

. 00 



O X 

o 



m^OQCQOQCQOQmOQ 



- 52 - 



- 53 - 



Table 4. Integrated Masses 



Conversion Factor^ log Si,^ log Mhi 

Galaxy (M0(Jy beam"^ m s~^)~^) (Jy beam~^ m s~^) (M0) 



Im Galaxies 



CVnIdwA 


52.98 


5.95 


7.67 


DDO 43 


595.15 


5.45 


8.23 


DDO 46 


530.30 


5.55 


8.27 


DDO 47 


134.65 


6.46 


8.59 


DDO 50 


125.49 


6.75 


8.85 


DDO 52 


1404.03 


5.28 


8.43 


DDO 53 


169.40 


5.49 


7.72 


DDO 63 


152.10 


6.01 


8.19 


DDO 69 


9.66 


5.86 


6.84 


DDO 70 


4.35 


6.97 


7.61 


DDO 75 


16.28 


6.65 


7.86 


DDO 87 


592.62 


5.62 


8.39 


DDO 101 


331.75 


4.84 


7.36 


DDO 126 


293.50 


5.69 


8.16 


DDO 133 


43.13 


6.39 


8.02 


DDO 154 


129.01 


6.35 


8.46 


DDO 155 


19.85 


5.70 


7.00 


DDO 165 


131.87 


6.01 


8.13 


DDO 167 


215.81 


4.84 


7.17 


DDO 168 


191.09 


6.17 


8.47 


DDO 187 


66.04 


5.30 


7.12 


DDO 210 


3.79 


5.72 


6.30 


DDO 216 


12.41 


5.65 


6.75 


F564V3 


351.11 


5.07 


7.61 


IC 10 


7.08 


6.93 


7.78 


IC 1613 


4.61 


6.86 


7.53 


LGS3 


2.09 


4.87 


5.19 


MSldwA 


124.28 


5.09 


7.18 


NGC 1569 


179.04 


6.14 


8.39 
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Table 4 — Continued 





Conversion Factor*^ 




log 


log Mhi 


Galaxy 


(M0(Jy beam-^ m s"^)~^) 


(Jy beam ^ m s ^) 


(Mq) 


NGC 2366 


132.41 




6.72 


8.84 


NGC 3738 


326.55 




5.54 


8.06 


NGC 4163 


69.27 




5.32 


7.16 


NGC 4214 


87.28 




6.82 


8.76 


Sag DIG 


6.47 




6.13 


6.94 


UGC 8508 


109.21 




5.24 


7.28 


WLM 


12.21 




6.76 


7.85 


BCD Galaxies 


Haro 29 


416.32 




5.18 


7.80 


Haro 36 


1010.34 




5.15 


8.16 


Mrk 178 


210.35 




4.67 


6.99 


VIIZw 403 


125.27 




5.59 


7.69 



^Factor to convert flux in the ROBUST-weighted integrated H I map to mass 
in solar masses. Note this conversion factor contains the necessary division 
by the number of pixels per beam, as discussed in §4. 

'^Integrated flux measured on the integrated Hi map made from the 
ROBUST-weighted image cube. 
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Table 6. Surface Density Photometry Parameters^ 



R.A. (2000.0) Decl. (2000.0) P.A.'' Incl 
Galaxy (hh mm ss.s) (dd mm ss) (deg) h/a'^ (deg) 



Im Galaxies 





1 9 '^8 4n 9 


4-'^9 d'^ 4n 


8n 

ou 


n 


78 

. 1 o 


41 n 

4:1 .U 




7 98 17 8 

/ ZrO -L I .O 


-l-4n 4fi 1 ^ 


u.o 




u. 


70 
. / u 


"4:0.0 


nun a(\ 


7 41 9fi fi 
/ ^± zu.u 


4-4n nfi 


84 
o^ 


n 

u. 


8Q 
.oy 




DUO A7 


7 41 i^i^l ^ 


4-1 fi 48 08 
T^±u uo 


70 
- / u 


n 

u. 


l^1 






R 1 Q DR 7 


-l-7n A'i 9^ 


1 8 

lO 


n 


79 


A(\ 7 






4-41 f^l 91 


K 
O 


n 


fi7 






8 ^4 08 n 


4-fifi in '^7 

n^uu lU Ol 


81 

Ol 


n 
u, 


1^1 




DDO fi^ 

VJ \-f\J uo 


Q 40 4 


4-71 1 1 09 

1 -L ±1 UZ 





1 

± . 


00 

.uu 



u.u 




Q i^Q 9^ n 


4-^n 44 49 


-u^ 




u. 


. ou 


uu.o 




1 n Dfi fin Q 


4- 1 Q "^n 

O -Lc/ fJU 


88 

OO 





5Q 
.oy 


^7 R 


DDO 7"^ 


in in "iQ 9 


n4 41 '^F, 

Ut: t:X OyJ 


41 


n 

u. 


8"^ 
.oo 


oo. (J 


DDO ^7 


1 n 4Q "^4 7 


QT 4fi 

n^UtJ ox 


7fi 


n 


.oo 


oo.u 


JJJJU iUi 


ii 00 oy.4 


1 Q1 Q1 nc 


-by 


n 
U, 


.by 


/I n /I 


DDO 126 


12 27 06.5 


+37 08 23 


-41 


0, 


.47 


67.7 


DDO 133 


12 32 55.4 


+31 32 14 


-6 


0, 


.69 


49.4 


DDO 154 


12 54 06.2 


+27 09 02 


46 


0, 


.50 


65.2 


DDO 155 


12 58 39.8 


+14 13 10 


51 


0, 


.71 


47.6 


DDO 165 


13 06 25.3 


+67 42 25 


89 


0, 


.54 


61.9 


DDO 167 


13 13 22.9 


+46 19 11 


-23 


0, 


.65 


52.8 


DDO 168 


13 14 27.2 


+45 55 46 


-24.5 


0, 


.63 


54.5 


DDO 187 


14 15 56.7 


+23 03 19 


37 


0, 


.80 


39.0 


DDO 210 


20 46 52.0 


-12 50 51 


-85 


0, 


.48 


66.9 


DDO 216 


23 28 35.0 


+14 44 30 


-58 


0, 


,45 


69.4 


F564-V3 


9 02 53.9 


+20 04 29 


8.5 


0, 


.83 


35.8 


IC 10 


20 21.9 


+59 17 39 


-38 


0, 


.78 


41.0 


IC 1613 


1 04 49.2 


+2 07 48 


71 


0, 


.81 


37.9 


LGS 3 


1 03 55.2 


+21 52 39 


-3.5 


0, 


.51 


64.4 


M81dwA 


8 23 57.2 


+71 01 51 


86 


0, 


.73 


45.8 


NGC 1569 


4 30 49.8 


+64 50 51 


-59 


0, 


.55 


61.1 
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Table 6 — Continued 





R.A. (2000.0) 


Decl. (2000.0) 


P.A.i' 




Incl 


Galaxy 


(hh mm ss.s) 


(dd mm ss) 


(deg) 


b/a" 


(deg) 


NGC 2366 


7 28 48.8 


+69 12 22 


32.5 


0.42 


72.1 


NGC 3738 


11 35 49.0 


+54 31 23 





1.00 


0.0 


NGC 4163 


12 12 09.2 


+36 10 13 


18 


0.64 


53.7 


NGC 4214 


12 15 39.2 


+36 19 38 


16 


0.91 


25.8 


SagDIG 


19 30 00.6 


-17 40 56 


87.5 


0.53 


62.7 


UGC 8508 


13 30 44.9 


+54 54 29 


-60 


0.54 


61.9 


WLM 


01 59.2 


-15 27 41 


-2 


0.44 


70.3 


BCD Galaxies 


Haro 29 


12 26 16.7 


+48 29 38 


86 


0.58 


58.6 


Haro 36 


12 46 56.3 


+51 36 48 


2 


0.81 


37.9 


Mrk 178 


11 33 29.0 


+49 14 24 


-50 


0.46 


68.6 


VIIZw 403 


11 27 58.2 


+78 59 39 


-10 


0.49 


66.0 



'^Prom the 1^-band (Hunter &; Elmegreen 2006). 
''Position Angle (P. A.). 

'^Minor-to-major axis ratio as measured in an outer F-band surface 
brightness isophote. For computing the inclination of the galaxy, 
the intrinsic h/a is assumed to be 0.3 for dwarf galaxies (Hodge & 
Hitchcock 1966, van den Bergh 1988). 



